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ABSTRACT 


This  report  is  a  continuation  of  a  study  of  the  effects  of  internal 
heat  transfer  on  the  temperature  of  hollow  spacecraft  and  the  require¬ 
ments  for  thermal  modeling.  Considered  herein  is  the  effect  of  in¬ 
ternal  heat  transfer  by  radiation  on  the  temperature  distribution.  The 
equation  governing  the  heat  transfer  of  a  spherical  shell  exposed  to 
parallel  radiation  is  derived;  conduction  and  radiation  are  considered. 
The  general  equation  is  simplified  by  assuming  steady  state,  and  a 
numerical  method  is  given  to  solve  the  steady  state  equation.  A  com¬ 
puter  program  is  described  which  employs  the  method.  Solutions  of 
the  steady  state  equation  are  graphically  presented  and  discussed. 

The  requirements  for  temperature  preservation  in  thermal  modeling 
are  derived.  The  possibility  of  thermal  modeling  without  temperature 
preservation  is  discussed.  It  is  observed  that  for  an  inside  emissivity 
to  outside  emissivity  ratio  greater  than  one,  the  requirement  for  dupli¬ 
cation  of  the  other  dimensionless  ratio  can  be  relaxed. 
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NOMENCLATURE 

Conducting  area 
Projected  area 
Radiating  area 
Radius  of  sphere 
Shell  thickness 
Specific  heat  capacity 
Density 

Ratio  of  emissivities 

Parallel  radiation  function 

Form  factor  from  area  at  9  to  area  at  9' 

Parallel  radiation  function 

Form-surface  factor  from  area  at  0  to  area  at  0' 
Step  function  to  limit  F{x) 

Step  function  to  limit  f (0) 

Thermal  conductivity 
Dimensionless  ratio 
Dimensionless  ratio 
Dimensionless  ratio 

Rate  of  heat  transfer  into  sector  by  conduction 

Rate  of  heat  transfer  out  of  sector  by  conduction 

Rate  of  heat  transfer  into  sector  from  source 

Rate  of  heat  transfer  out  of  sector  by  radiation 

Rate  of  heat  transfer  into  sector  by  radiation 

Net  rate  of  heat  transfer  into  sector 

Parallel  radiation  function 

Area  average  of  pc 

Distance  from  axis  of  symmetry 

Temperature 

Arbitrary  temperature 
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CT 


Time 

Arbitrary  time 
Volume 
Cosine  of  6 
Step  size 

Dimensionless  time 
Dimensionless  temperature 
Lower  bound  of  Z(-l) 

Trial  value  of  Z(-l) 

Upper  bound  of  Z(-l) 

Inside  absorptivity 

Absorptivity  of  radiation  from  source 

Inside  emissivity 

Outside  emissivity 

Angle  from  axis  of  symmetry 

Inside  reflectivity 

Stefan- Boltzmann  constant 
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•  SECTION  I 
INTRODUCTION 


This  report  is  a  continuation  of  a  study  of  the  effects  of  internal 
heat  transfer  on  the  temperature  of  hollow  spacecraft  and  the  require¬ 
ments  for  thermal  modeling.  Considered  herein  is  the  effect  of  inter¬ 
nal  radiative  heat  transfer  on  the  temperature  distribution. 

In  Ref.  1,  the  effect  of  internal  convection  on  the  temperature  of 
a  spacecraft  model  of  arbitrary  shape,  subjected  to  parallel  radiation, 
was  considered,  and  the  transient  temperatures  were  calculated.  With 
the  aid  of  numerical  results,  the  conditions  under  which  convection  can 
be  neglected  were  determined.  Also,  thermal  modeling  rules  were 
derived  for  testing  scale  models.  Thermal  modeling  is  a  valuable  tech¬ 
nique  in  ground  testing  of  spacecraft.  Several  aspects  of  thermal  model 
testing  in  space  simulation  chambers  are  discussed  in  Ref.  2. 

In  this  report,  the  geometry  is  restricted  to  a  spherical  shell,  and 
solutions  are  obtained  for  the  steady  state  case.  The  spherical  shell  is 
subjected  to  parallel  radiation  which  is  considered  uniform  for  the  cal¬ 
culations  presented  herein;  however,  the  equations  derived  and  the 
numerical  method  allow  axially  symmetric  nonuniformities.  These 
somewhat  arbitrary  restrictions  were  imposed  so  that  relatively  simple 
calculations  resulted,  which  nevertheless  bring  out  features  of  a  more 
general  nature.  By  carrying  out  selected  calculations  of  this  type,  it 
is  hoped  that  better  insight  can  be  obtained  concerning  the  importance 
of  various  parameters  to  thermal  testing  in  space  simulation  chambers. 


SECTION  II 

MATHEMATICAL  ANALYSES 


2.1  GENERAL  EQUATIONS 

The  system  to  be  considered  is  shown  in  Fig.  1  with  vacuum  inside 
and  outside  of  the  sphere.  It  is  assumed  that  the  shell  thickness,  b, 
is  small  enough  to  disallow  any  temperature  gradient  in  the  radial  direc¬ 
tion.  It  is  also  assumed  that  the  inside  of  the  sphere  emits  and  reflects 
diffusely.  The  parallel  radiation,  ps,  is  allowed  to  be  a  function  of  r. 
This  may  be  useful  at  a  later  time  in  investigating  the  effects  of  non¬ 
uniformities  of  solar  simulators. 
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The  general  equation  governing  the  system  can  be  obtained  by  per¬ 
forming  a  heat  balance  on  the  sector  of  the  sphere  between  6  and  6  +  d£> 
and  the  equation 


dpi„  =  cDdvfl 
a  t 


(1) 


The  rate  of  heat  transfer  into  the  sector  at  6  by  conduction  (Fig.  1) 


is 


Pi 


=  — kAk  -1  *1 

a  86 


The  rate  of  heat  transfer  out  of  the  sector  by  conduction  at  0  +  d 0  is 


■  a  k  1  0 

*  -kA  ;  m 


s+iie 


=  _kAk  I  d-l  _  Js_i_ 

a  dd  a  89 

•$+de 


(A‘  ID" 


Thus  the  total  rate  of  heat  transfer  into  the  sector  by  conduction  is 


(2) 


Let  q  be  the  average  over  the  projected  area  of  the  sphere  of 


p  •  then 


or 


qs 


qs 


=  Tpf  Ps(r)  2,7  r 

**  o 

=  ~T  I  r  ps(r)  dr 

a  •'o 


dr 


(3) 


Now  define  the  dimensionless  functions 

=  -jp  ps  (a  sin  B) 


and 


i1  o  <  0  <  f 

h(0)  =  < 

/o  f  <  0  <  - 


(4) 


(5) 


The  rate  of  heat  transfer  into  the  sector  from  the  parallel  radiation 
source  is  then 

dpi  =  as  ps(r)  dAp 


or 


dp3  =  asqs  f(0)  h(0)  dA1* 


(6) 
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The  rate  of  heat  transfer  emitted  out  of  the  sector  is  given  by 

dp4  =  (<j  -  <0)o  T4dAR  (7) 

The  rate  of  heat  transfer  into  the  sector  from  radiation  inside  the 
sphere  is  given  by 

dp,  =  ai  JntiovmdGe'..oAAe'  (8) 

0 

where  the  integration  is  of  the  area  at  6'. 

The  different  areas  are  given  by 

A*1  =  2rr  a  b  sin  9  (9) 

dA^  =  tt  a2  d  (sin2  9) 
or 

dAP  =  2  it  a2sin0cos#d0  (10) 


The  volume  is 


dAR  s=  2 rr  a2  sin  9  d9 


dV  =  b  dA* 


dV  =  2/r  a2  b  sin  0  d0 


It  is  proved  in  Appendix  I  that 

dG6'-9  =  sin  0 

A  substitution  in  Eq.  (1)  of  dp^n  =  (pj  -  pg)  +  dp3  -  dp4  +  dp5 
results  in 

2,rab±-[cos0g  +  Sin0 


+  2jra2osqs  f(0)  h(0)  sin  9  cos  eae 


-  2na2  (ci  ■*.  e0)  a T4  sin  Odd 


a;  f;  a  sin  9  dfl^irn2  J  T+(0')  sin  d'dd' 


=  2j7S2  b  c  D  sin  9  d9 

a  t 
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Dividing  by  -Tsin0d0  for  0 this  becomes 


H1  6  M  ‘  S*] h  2d*  “-‘J-  f(6)  h(0)  cos  6 

/* 

T\d')  sin  0'  d0' 


=  -Ja2  b  c  I) 

01 


Define  the  dimensionless  variables 


Z  =  T 

*  m 


where,  for  the  present,  Tm  and  tm  are,  respectively,  some  arbitrary 
temperature  and  some  arbitrary  time.  Equation  (14)  is  obtained  in 
dimensionless  form  by  using  Eqs.  (15)  and  (16)  and  dividing  through  by 
2bkTm.  The  result  is 


coi  ed-l  +  ^  + 

do  de 2 


\bkTm/ 


f(0)  h(0)  cos  0 


-  ('  *  r)2‘ +  fe)/z,(n,to rir  (17) 

O 

_  /a2  c  P\  dZ 
\  k  ln  /  <?v 

This  introduces  three  dimensionless  quantities  which  will  be  defined 

Ns  =  /i  tn 

bkTm  (18> 


Nn  = 


a2f0<^Tm 


These  quantities  may  be  used  as  a  measure  of  the  effects  of  an 
external  energy  source,  external  radiation,  and  the  heat  capacity, 
respectively,  as  compared  with  the  heat  transfer  by  conduction.  Also 
define 

E  =  t  (2D 
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In  these  terms,  Eq.  (17)  becomes 


a2/ 

aa2  f 


C  Ol  V  -rr: 

dO 


Ns  m  MO)  cos  6 


+ 


Nr 


')  sin  d'<\Q' 


(L  + 


(22) 


=  N 


c 


az 

dy 


This  equation  can  be  written  in  a  different  form  by  a  change  of 
variable 


X  =  cos  6 


From  the  equation 


is  obtained 


Also 


Defining 


and 


a?  =  dz  Ax  =  _sin0  az 
aa  dx  de  dx 


cot  9  2Z  =  -  cos  a  p  =  -  X  & 


aa 


dx 


ilz_=  ±(.  .inert) 

dx2  aa\  dx) 


=  -  cos  a  —  +  sin2  a 
dx 


d'z 


=  -  x  S*  +  (1  -  x2) 

dx  dx 2 


F(x)  =  f  (cos  1  xl 


ll  0  <  x  <  l 


H(x)  =< 


0  -l  <  x  <  0 


Thus  Eq.  (22)  becomes 


(1  -  x2)  £%  -  2x  ^  +  Ns  F(x)  H(x)  x 
dx  dx 


d7. 


+  Nil 


[lEf  ZTxOdx'  -  (1  +  E) 


(24) 


(25) 


(26) 


Either  Eq.  (22)  or  Eq.  (26)  describes  the  temperature  of  the  system 
as  a  function  of  position  and  time.  Given  an  initial  temperature  distribu¬ 
tion,  the  numerical  solution  of  these  equations  could  be  obtained. 
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2.2  STEADY  STATE  EQUATIONS 

If  steady  state  condition  is  assumed,  a  simplification  of  Eqs.  (22) 
and  (26)  can  be  made.  Performing  a  heat  balance  on  the  whole  sphere 
at  steady  state  gives 


e0  a  T4  dAR  -  ne2  as  q« 


>  a  J*  T4  sin  6  A9  =  na 1  as  qs 


from  which  is  obtained 


sin  0  dO  =  Qs  .^a 
2en  a 


f  Z4  sin  6  <16  =  — °s  (29> 

**o  2(oa  I  m 

If  Tm  is  defined  as  the  temperature  such  that  if  the  sphere  was  isother¬ 
mal,  it  would  radiate  the  same  power;  then 


/” 

e„  a  T4  dA« 


T  - 

A  ni  — 


With  this  definition,  Eq.  (29)  becomes 


Thus 


Z4  sin  0  Ad  =  2 


Z4dx  =  2 


Also,  for  steady  state  conditions  a  relation  between  Ng  and  Nr  can  be 
obtained.  By  Eq.  (18) 


v  _  »afl3  qB  T^, 

■  s  b  k  TjJ, 

qs  Ta 

bKplM 

\4^0-/ 


4aJ  ffT; 
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or 

Ns  =  4Nr  (33) 

Using  these  substitutions  and  noting  that  for  steady  state  the  right 
hand  side  is  zero  and  the  partial  derivatives  become  total  derivatives, 
Eq.  (22)  becomes 

0  +  cot  6  ||  +  4N|j  f(0)  h(0)  cos  6 

+  XRg)r,  (2)  -  Nr  (l  +  E)  Z4  =  0 
or 

+  cot  e  &  +  nr  |4f(e)  h(e)  cos  e+  e  -  a  +  e)  z4j  =  o  (34) 
Similarly  Eq.  (26)  becomes 

(1  -  xJ)  -  2x|-z+  Nr  |4F(x)  IF(x)  x  t  E  -  (1  +  E)Z41  =  0  <j3 5) 

l  1 

Along  with  the  simplification  in  the  equations  is  a  complication  in 
the  condition  to  be  satisfied  by  the  solution.  Whereas  the  condition  to 
be  satisfied  by  the  solution  of  the  transient  equations  was  simply  the 
initial  temperature  distribution,  the  solution  of  Eq.  (34)  must  satisfy 
Eq.  (31)  and  the  solution  for  Eq.  (35)  must  satisfy  Eq.  (32).  This  ■ 
complication  in  side  conditions  is  not  peculiar  to  this  particular  sys¬ 
tem  but  is  a  general  circumstance  when  going  from  the  transient  equa¬ 
tions  to  the  steady  state  equations  of  most  systems. 


SECTION  III 

SOLUTION  OF  THE  STEADY  STATE  EQUATIONS 


3.1  NUMERICAL  METHOD 

This  section  describes  the  numerical  method  used  to  obtain  numeri¬ 
cal  solutions  to  the  steady  state  equations.  A  steady  state  solution  could 
be  obtained  by  assuming  an  initial  temperature  distribution  and  obtaining 
a  numerical  solution  of  the  transient  equations.  However,  a  more  effi¬ 
cient  method  is  to  assume  a  starting  value  of  the  steady  state  equation, 
obtain  a  solution,  and  iterate  to  find  the  solution  which  fits  the  side  condi¬ 
tion.  This  method  is  described  below  as  applied  to  Eq.  (35). 
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Define 

and  Eq.  (35)  becomes 

.is 


S  =  i* 
dx 


(36) 


(1  -  x2) -  2x S  +  Nn  V‘(x)]!(x)x  +  E  -  (L 

dx  ( 


+  E)  7.*j  =  0 


Also  define 


L* 


R(x)  =  /  Z4  d> 


From  these  equations,  we  easily  obtain  the  system  of  equations 


d  x 


dS  =  l 

dx  I  —  x2 


|  2x  S  +  N||  Jj(l  +  E)Z4  -  4F(x)  H(x)  x  -  eJ| 


d  y.  =  c 
dx  ™  1 


(37) 

(38) 

(39) 

(40) 


A  numerical  solution  of  this  system  of  equations  could  easily  be  obtained 
if  the  values  of  R,  S,  and  Z  at  x  =  -1  were  known.  However,  instead  the 
following  side  conditions  must  be  met.  From  Eq.  (37) 

R(-l)  =  0  (41) 

From  Eq.  (32)  is  obtained 

RU)  =  2  (42) 

and  from  Eq.  (35)  is  obtained  the  condition 


S(-l)  =  fNn  {(1  -  E)[z(-i)]4  -e} 


(43) 


From  the  physical  orientation,  the  coldest  spot  on  the  sphere  is 
at  x  =  -1;  thus.S(-l)  is  positive  or  zero  and  from  Eq.  (43)  is  obtained 

ZM)  *  (ttk) 

Also  with  this  assumption,  from  Eq.  (32), 

Z(-l)  <  1 

If  and  Zu  are  lower  and  upper  bounds  of  Z(-l)  then  by  the  above  in¬ 
equalities  valid  values  are 


=  (*) 
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At  the  beginning  of  each  iteration  a  value,  ZQ,  is  defined  as 

Z°  =  |(Zj  +  Zu) 

and  starting  values 

R(-l)  =  0 

S(-l)  =  {Nr  [(1  -  E)  Z40  -  E] 

Z(-l)  =  Z0 

i 

are  used.  It  is  seen  that  the  side  conditions,  Eqs.  (41)  and  (43),  are 
satisfied.  Iteration  is  continued  until  Eq.  (42)  is  satisfied  within  a 
given  tolerance,  AR.  That  is  until 

2  -  AR  <  R(l)  <  2  +  AR 


Each  iteration  determines  whether  Z0  is  an  upper  or  lower  bound 
for  Z(- 1).  It  is  assumed  that  if  ZQ  is  too  low,  then  the  solution  obtained 
for  Z(x)  is  below  the  correct  solution  and  that  if  ZQ  is  too  high,  then  the 
solution  obtained  is  above  the  correct  solution.  (This  assumption  was 
verified  by  numerical  results. )  If  at  any  point  in  a  solution 


ZU)  <  Z t 

then  it  is  known  that  ZQ  is  a  lower  bound  for  Z(-l).  Therefore  Zjg  is  set 
equal  to  ZQ  and  another  iteration  is  begun.  If  at  any  point  in  a  solution 

RU)  >  2  +  AR 

then  it  is  known  that  Z0  is  an  upper  bound,  so  Zu  is  set  equal  to  ZQ  and 
a  new  iteration  is  begun.  If  a  solution  proceeds  to  X  =  1,  then  R(l)  is 
tested  to  see  if 

R(l)  <  2  -  AR 

If  so,  then  Z^  is  set  equal  to  ZQ  and  another  iteration  begun.  If  not, 
then  Eq.  (42)  is  satisfied  within  the  given  tolerance  and  the  desired 
solution  has  been  obtained. 

A  flow  chart  of  this  method  is  shown  in  Fig.  2. 
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Fig.  2  Flow  Chart  of  Main  Program 
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3.2  COMPUTER  PROGRAM 

A  computer  program  was  written  in  FORTRAN  II  for  the  SDS  920 
computer  to  solve  the  steady  state  equation  employing  the  method 
described  above.  The  main  program  follows  the  flow  chart  shown  in 
Fig.  2.  The  subroutine  STEP  called  by  the  main  program  is  a  stand¬ 
ard  subroutine  used  to  solve  differential  equations.  STEP  calls  an¬ 
other  subroutine  FUN  which  evaluates  Eqs.  (38),  (39),  and  (40).  These 
are  explained  in  detail  below. 

3.2.1  Main  Program 


The  input  to  the  program  is  N,  MAX,  and  L  defined  below  and 
Nr,  E,  and  AR.  N  is  the  number  of  steps.  Since  the  solution  is 
required  for  -1  <  x'<  1,  the  step  size  is 

Ax  =  2/N 

MAX  contains  a  limit  to  the  number  of  iterations.  The  fractional  por¬ 
tion  of  the  SDS  920  floating  point  number  is  39  binary  digits.  Since  the 
numerical  method  halves  the  interval  containing  Z(-l)  each  iteration, 
the  computer  precision  limits  the  method  to  39  iterations.  If  L,  is  1, 
then  the  solution  is  printed  after  every  iteration.  If  L  is  2,  then  only  the 
last  solution  is  printed. 

The  flow  chart  for  the  main  program  is  shown  in  Fig.  2,  and  the 
listing  is  given  in  Appendix  II. 

3.2.2  Subroutine  STEP 


Subroutine  STEP  is  a  program  of  the  Runge-Kutta  one  step  method 
of  solving  differential  equations.  This  method  is  explained  in  Ref.  4. 


Given  a  system  of  differential  equations 


dy‘ 


=  fl(x,  y*.  y3 . yn),  i  =  1,  2, 

dx 


and  starting  values 


>rl(x0) 


>0 


(44) 


A  one  step  method  is  an  algorithm  which  uses  the  differential  equations 
and  the  starting  values  to  find  an  approximation  to  the  solution  at 
xQ  +  Ax 

y'(x0  +  Ax)  =  yl 
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These  values  can  be  used  as  new  starting  values  to  find  an  approxima¬ 
tion  to  the  solution  at  x0  +  2Ax  or  in  general,  after  p  steps 

y'(xp)  =  vp 

where 

Xp  =  x0  +  p  Ax 

The  smaller  the  step- size,  Ax,  is  the  better  the  approximation  be¬ 
comes.  The  accuracy  is  limited  only  by  the  round  off  error  which 
depends  on  the  precision  of  the  computer. 

Subroutine  STEP  is  called  by  the  statement 

CALL  STEP  (N,  X,  Y,  DX) 

where  Y  is  dimensioned  Y{25).  N  contains  the  number  of  equations  and 
DX  the  step  size.  Suppose  p  steps  have  been  taken.  Then,  when  STEP 
is  called,  X  contains  Xp  and  Y  contains  yp,  with  i  =  1,  2,  ...  n.  On 

return,  X  contains  Xp+j  and  Y  contains  y^  j,  with  i  =  2,  ...  n.  STEP 
calls  on  a  subroutine  FUN  to  evaluate  Eq.  (44). 

The  listing  of  STEP  is  given  in  Appendix  II. 


3.2.3  Subroutine  FUN 

In  the  present  case,  subroutine  FUN  was  written  to  evaluate  Eqs. 
(38),  (39),  and  (40)  which  correspond  to  Eq.  (44)  for  this  problem.  It 
is  seen  that  Eq.  (39)  is  singular  at  x  =  ±1.  To  avoid  this  difficulty 
ifx<-l  +  0.1  Ax,  then  it  is  assumed  that 

1  -  xJ  =  1  -  x\ 

where  xj  =  -1  +  0.  1  Ax.  If  x  >  1  -  0.  lAx,  then  the  above  approximation 
is  again  used  only  this  time  xj  =  1  -  0.  lAx. 

The  flow  chart  of  FUN  is  shown  in  Fig.  3,  and  the  listing  is  given 
in  Appendix  II. 


SECTION  IV 

RESULTS  AND  DISCUSSION 


The  computer  program  described  above  was  used  to  obtain  solu¬ 
tions  of  Eq.  (35)  for  various  values  of  N^  and  E.  The  convergence 


13 


AEDC-TR-67-254 


Fig.  3  Flow  Chart  of  FUN 
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of  the  iterations  is  illustrated  in  Fig.  4.  Figures  5,  6,  and  7  are 
the  graphs  of  solutions  with  Nr  =  1,  2,  and  3,  respectively,  for 
various  values  of  E.  The  solution  with  Nr  infinite  can  be  obtained 
from  Eq.  (35),  where  for  this  case,  7A  is  a  linear  function  of  x.  This 
linearity  is  shown  in  Fig.  8,  and  the  Z  versus  x  relationship  is  shown 
in  Fig.  9.  The  case.  Nr  =  0,  can  arise  only  from  infinite  conductivity 
since  it  is  assumed  that  b  «  a.  This  means  the  sphere  would  be  iso¬ 
thermal  giving  the  solution,  Z(x)  =  1,  independent  of  E.  It  is  interest¬ 
ing  to  note  that  at  x  =  0.  25,  the  dimensionless  temperature  Z  is  very 
close  to  one  in  all  cases.  The  solutions  show  that  when  E  >  1,  the 
temperature  distribution  changes  relatively  slow  as  Nr  is  changed. 
That  is,  the  temperature  distribution  is  determined  largely  by  internal 
radiation.  For  E  <  1,  Nr  has  a  greater  influence  on  the  solutions, 
implying  that  conduction  also  is  an  important  factor  in  determining  the 
temperature  distribution. 


If  Nr  and  E  are  the  same  for  two  different  systems,  then  the  same 
equation  describes  both  systems;  thus,  the  temperature  distribution  of 
one  system  can  be  inferred  from  the  measured  temperature  distribution 
of  the  other  system.  This  is  the  basis  for  thermal  modeling.  The  ob¬ 
servations  made  above,  on  the  behavior  of  the  solutions  of  Eq.  (35), 
imply  that  for  E  >  1  the  tolerance  of  duplication  of  Nr  in  thermal  model¬ 
ing  can  be  relaxed,  but  when  E  <  1  the  requirement  for  duplication  of 
Nr  becomes  more  stringent. 


In  thermal  modeling  it  may  be  desired  to  preserve  temperature, 
since  the  thermal  properties  may  be  a  function  of  temperature.  For 


this  case,  in  addition  to  preserving  Nr  and  E,  one  must  preserve  Tm 
which  is  given  by  Eq.  (30).  Thus,  from  Eqs.  (19),  (21),  and  (30), 
after  cancelling  constants  and  Tm,  are  obtained  the  requirements 

&).  ■ 

(45) 

(t).  ■  a 

(46) 

s(la\  /®  s  (1  s\ 

\  f“',.  ~  V  f”  />> 

(47) 

where  the  subscripts  m  and  P  refer  to  model  and  prototype.  If  further¬ 
more,  materials  and  surface  properties  are  preserved,  then  Eqs.  (46) 
and  (47)  are  automatically  fulfilled  and  from  Eq.  (45)  the  well-known 
scaling  requirement 


is  obtained. 


(f)  - 


(48) 
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1  E  -  0  AR  =  0.01 


N  =  40  NR  - 
12th  Iteration  Was  within  Tolerance  - 


X 


Fig.  5  Solutions  with  Nr  =  1 
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Fig.  6.  Solutions  with  NR 
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Fig.  7  Solutions  with  Nr  «•  3 
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There  is,  at  least  theoretically,  a  possibility  of  modeling  without 
temperature  preservation.  The  preservation  of  Np  requires  that 


But  by  Eq.  (30) 

(T-)- 
(T">  ' 


(49) 


(50) 


Thus,  from  Eqs.  (49)  and  (50)  is  obtained  the  requirement 


s' 
3 _ 1 

qs 

,  \ 
q&  *o  'p 

The  preservation  of  E  requires  that 


(51) 


(52) 

(fi)p  P 

Given  a  prototype,  after  a  suitable  choice  of  dimensions  and  selection 
of  material  for  the  model,  the  left-hand  side  of  Eq.  (51)  is  determined. 
The  choice  of  surface  finish  on  the  model  will  determine  (Oglm  and 
(e0)m.‘  It  is  then  possible  to  adjust  (qs)m  to  fulfill  Eq.  (51).  Then 
(e^m  is  determined  by  Eq.  (52).  Once  the  temperature  distribution  of 
the  model  is  measured,  the  temperature  of  the  prototype  can  be  calcu¬ 
lated  by  Eq.  (15). 


Wp  (TL 


(53) 


since  Z  is  preserved  by  design.  Equation  (30)  is  used  to  calculate  Tm. 


SECTION  V 
SUMMARY 


The  equation  governing  the  heat  transfer  of  a  spherical  shell  sub¬ 
jected  to  parallel  radiation  was  derived,  conduction  and  radiation  being 
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considered.  A  numerical  method  was  developed  to  solve  the  steady  state 
equation,  and  a  computer  program  was  described  which  employed  the 
method.  Solutions  of  the  steady  state  equation  were  graphically  presented. 

The  requirements  for  temperature  preservation  in  thermal  model¬ 
ing  were  derived.  The  possibility  of  thermal  modeling  without  tempera¬ 
ture  preservation  was  discussed.  It  was  observed  that  for  an  inside 
emissivity  to  outside  emissivity  ratio  greater  than  one,  the  require¬ 
ment  for  duplication  of  the  other  dimensionless  ratio  can  be  relaxed. 
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APPENDIX  I 

DERIVATION  OF  FORM-SURFACE  FACTOR 


The  fact  that  flux  leaving  any  part  of  the  interior  of  a  sphere  diffusely 
spreads  uniformly  over  the  interior  of  the  sphere  implies  that  the  form 
factor  from  the  area  at  0  to  the  area  at  0' 


F6-6' 


dAj- 

4ffaJ 


Also,  the  form  factor  from  the  sphere  to  the  area  at  0'  is 


Fs-fl'  = 


dA  g ' 
in  a1 


Consider  flux  leaving  the  area  at  0.  ^q.q>  of  this  flux  reaches  0' 

directly.  The  flux  from  any  reflection  leaves  the  sphere  uniformly  so 
Fg _q,  of  the  reflected  flux  reaches  0'. 


Of  the  flux  leaving  0,  the  flux  eventually  striking  0'  is,  by  adding 


G 0_0'  =  Fe-d'  +  Pi  Fs-0'  +  p-  Fs-0'  + 


dA 


477  a 


>( 


1  *•  Pi  +  Pi 


>  +  *  * 


R 

dA  6' 


477 aJ  (1  -  pj) 


or  from  Eq.  (11)  this  becomes 

Gd_d>  =  h£C 

This  is  the  fraction  of  the  flux  leaving  the  area  at  0  which  eventually 
strikes  the  area  at  0'  and  is  called  the  form- surface  factor. 
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APPENDIX  II 

FORTRAN  LISTING  OF  COMPUTER  PROGRAM 


■ 

1 

C 

5027  RADIATING  SPHERE  D.C.  TB DD  1-12-67 

■ 

2 

DIMENSION  XLlOOJjZItQQl.ZFl  100I.RI  1G0I.SI  lOQi.DSllOOi.Vl3  l.DUM 

■ 

3 

COMMON  D.E.DX.ZFO.DSO 

■ 

4 

_C 

■ 

5 

C 

INPUT  AND  INITIATE 

■ 

6 

C 

■ 

7 

1  ACCEPT  TAPE  1001.N.MAX.L 

■ 

9 

IFINI2.2.3 

■ 

g 

2  PRINT  1004 

■ 

10 

STOP 

■ 

u 

3  ACCEPT  TAPE  1002.D.E.T 

■ 

_ 12 

DX-2./N 

■ 

13 

XII  J--1. 

■ 

14 

ruj-oa 

■ 

15 

ZU-1. 

■ 

16 

ZL-IE/I 1,  +  El I**, 25 

■ 

17 

ISUIT-2 

■ 

19 

IT-o 

■ 

19 

C 

■ 

20 

c 

START  OF  ITERATION 

■ 

21 

c 

■ 

22 

4  I  ?  - 1  T  M 

■ 

23 

Z0-.5M  ZU+ZLI 

■ 

24 

Sill -.5* DM C l.*El*Z0**4-El 

■ 

25 

Zi  1 1-ZO 

■ 

26 

U--1. 

■ 

27 

vm-o. 

■ 

26 

VI2I-SI 1  I 

■ 

29 

VI3I-Z0 

■ 

30 

CAUL  FUNIU.V.DUM) 

■ 

31 

ZFI  1  I-ZFO 

■ 

32 

DSllJ-DSO 

■ 

33 

c 

FIND  SOLUTION 

■ 

34 

DO  5  I-l.N 

■ 

35 

CALL  STEP l 3.U. V, DX 1 

■ 

36 

XI  1*1 I-U 

■ 

37 

ZI  IMI-VI3I 

■ 

36 

ZFI I'M  I-ZFO 

■ 

39 

R 1  I  ♦  1 1  ■  V  I  1  I 

■ 

40 

SI  U1I-VI2I 

■ 

41 

DSI 1*1 I-DSO 

■ 

42 

IFIVI3I-ZU50.50.60 

■ 

43 

60  C-2.-V11I 

■ 

44 

|FIT*C 151.51.5 

• 

45 

5  CONTINUE 

■ 

46 

c 

■ 

47 

c 

TEST  TO  SEE  IF  THIS  IS  LAST  ITERATION 

■ 

46 

c 

■ 

49 

I F ( T-C 1 50.6.6 

■ 

50 

6  IGUIT-3 

■ 

51 

SO  TO  10 

■ 

52 

50  ZL-ZO 

■ 

53 

GO  TO  52 

■ 

54 

51  ZU-ZO 
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■ 

55 

52 

Z 1  I *2  J “0 

■ 

56 

7 

IF1MAX-IT18.8.9 

■ 

52 

A 

ICUI T-l 

■ 

56 

CO  TO  10 

■ 

59 

9 

CO  T0(10#4|,L 

■ 

60 

C 

• 

61 

C  ST  ART  OUTPUT 

■ 

62 

C 

■ 

61 

_ 1_Q_ 

PRIM  1000 

■ 

64 

PRINT  1005 

■ 

65 

PRINT  1003,D«E,IIX.L»2L.ZU 

■ 

66 

PRINT  1006.IT 

■ 

67 

CO  T0L11, 11,121, I  QUIT 

■ 

66 

11 

PRINT  1008 

■ 

69 

CO  TO  13 

■ 

70 

12 

PRINT  1007 

■ 

71 

13 

PRINT  1009 

■ 

72 

PRINT  1003. IXt I ]. I-1.N*1 1 

■ 

73 

PRINT  1010 

m 

74 

PRINT  1003. IZtll. 1-1. N*ll 

m 

_ 7_5_ 

PRINT  1011 

m 

76 

PRINT  1003.IZF! I l.I-l.N*! ) 

m 

77 

PRINT  1012 

m 

76 

PRINT  1003. IRC  I ]. I*1.N»1 I 

m 

79 

PRINT  1013 

• 

60 

PRINT  1003.  IS!  I  1.  I-l.NM  J 

m 

61 

PRINT  1014 

m 

62 

PRINT  1003. I D S I  I 1  »  I  ■ 1 »  N»1 ) 

m 

83 

CO  mi>4,U,L5UiT 

m 

64 

1000 

F6RMATIS1D1C.  TODD  5027*1 

m 

85 

1001 

FORMA  T I 40 1 2 1 

m 

66 

1002 

FORM ATI 6E 12.01 

• 

6  7 

HI  I'M 

FORMAT! 1P10E12. 41 

m 

66 

1004 

FORMAT! 1 H 1 1 

m 

69 

1005 

FORMAT  I JH04X2HNR10X1HE11X2HDX10X1HTUX2H2L10X2HZU/1 

m 

90 

1006 

FORMAT! SOI TERAT IONS  131 

m 

91 

1007 

F0RM4TIS0*]THIN  TOLERANCES) 

m 

92 

1008 

FORMATISONOT  WITHIN  TOLERANCES! 

m 

_ 93 

mi  i*n 

FORMAT ( SOXS/  1 

m 

94 

1010 

FORMAT ! SOZS/ J 

m 

95 

1011 

FORMAT  ISOFOURTH  POWERS/1 

m 

96 

1012 

FORMAT  I SOINTECRAL  OF  FOURTH  POWERS/) 

m 

97 

1013 

FORMATISOFIRST  DERIVATIVES/! 

m 

96 

1014 

FORMAT  I SOSECOND  DERIVATIVES/! 

80HHBB9ALL0CATIBND 


77776  D 

77774  E 

77772  DX 

77770 

ZFO 

77766  DSO 

PROGRAM  ALLOCATION 

00004  X 

00314  Z 

00624  ZF 

C 1 1 34 

R 

01444  S 

01754  DS 

02264  V 

Q2272 

DUM 

02300  N 

02301  MAX 

02302  L 

02303 

IOUIT 

Q2304  IT 

02305  L 

02306  T 

02310 

ZU 

02^12  ZL 

02314  Z£ 

02316  U 

02320 

C 

SUBPROGRAMS 

REgUiRED 

FUN 

THE  FND 

STEP 
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1  C  ONE  STEP  8 Y  RUNGE-KUTTA  D.C.  TODD  1- 


■  2  SUBROUTINE  STEP  1 N. X*  Y  »  DX  J 

■  3  DIMENSION  Y1251,Dl25l.Y0(25].FI25l 

■  4  CAUL  FUMX.Y.DI 

-  5  DO  1  1-1. N 

«  6  YOm-YlJl 

-  7  Fill-Din 

_  a  i  yi  n-yom*.5*Dx*DC n 

■  9  X-X*.5*DX 

-  10  CALL  FUMX.Y.DI 

■  11  DO  2  1-1 .N 

-  12  fi  n-Fi  n*2.*Dm 

-  13  2  Yl) 1-YOI )  )♦  .5«DX*D I  1  I 

_ a _  14  _  CALL  FUMX.Y.DI 

-  15  DO  3  1-1. N 

-  16  FI)1-Fm+2.*D()) 

-  1 7  3  Yl 11-YOI  I  1  *DX *D (  I  1 

-  18  X-X+*5*DX 

-  19  CALL  FUMX.Y.DI 

■  _ 20  DO  4  I-l.N 

-  21  F| | 1-,166666666667-lFI |  I  ♦ D I ( (  1 

-  22  4  Y| 1  1-YOI  1  1  ♦ n X * F (  1  ] 

-  23  RETURN 

-  24  END 

PROGRAM  ALLOCATION 

Dummy  Y  00017  D  00101  YO 

00345  1  DUMMY  N  00246  STEP 

DUMMY  DX 


SUBPROGRAMS  REQUIRED 
FUN 

THE  FND 


12-67 


00163  F 
DUMMY  X 
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■ 

l 

C  DERIVATIVES  FOR  THE  RADIATInG  EQUATIONS  D.C.  TODD  1-12-67 

■ 

2 

SUBROUTINE  FUNtX.Y.Vl 

■ 

3 

DIMENSION  YI3J.V131 

■ 

4 

COMMON  JJ»£a1X,ZF.J1S 

■ 

5 

Xi»-1 .  ■*  «  1  *DX 

■ 

6 

IFIX1-XI1*3.3 

■ 

7 

1 

IF { X  12*2*4 

■ 

s 

_ 2_ 

Xt-X 

■ 

9 

3 

H-0. 

■ 

10 

GO  10  6 

■ 

11 

4 

H»1 . 

■ 

12 

Xi»l.-.1*DX 

■ 

13 

IFIX1-X 16*6*5 

■ 

14 

5 

Xl-X _ 

m 

1 5 

6 

2F- Y l 3 1 »  *4 

m 

16 

Vlll-ZF  _  .  . 

m 

1 7 

DS»12.«X*Y121*D«1 Il.+E1«ZF-4.*FIX1«H*X-El l/ll.-Xl*Xll 

m 

18 

y|2I-CS 

u 

19 

VI31-Y121 

u 

20 

RETURN 

■ 

21 

END 

COMMON  ALLOCAT ION 


77776 

D 

77774 

E 

77772 

DX 

77770 

IF 

77766 

DS 

PROGRAM 

ALLOCATION 

DUMMY 

Y 

DUMMY 

V 

00007 

FUN 

00011 

XI 

■>]iy;ii 

_X _ 

M'I'IU 

_ 

subprograms  REQUIRED 


F 


THE  END 


1  FUNCTION  F 1 X  1 

2  _ F»l. _ 

3  RETURN 

4  _ END _ 


PROGRAM  ALLOCATION 


00002  F 

The  end 
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